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Abstract 

We study the effective conductivity cTe for a random wire problem on the d- 
dimensional cubic lattice Z'^, d > 2 in the case when random conductivities on bonds 
are independent identically distributed random variables. We give exact expressions 
for the expansion of the effective conductivity in terms of the moments of the disorder 
parameter up to the 5th order. In the 2D case using the duality symmetry we also 
derive the 6th order expansion. We compare our results with the Bruggeman approx- 
imation and show that in the 2D case it coincides with the exact solution up to the 
terms of 4th order but deviates from it for the higher order terms. 
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1 Introduction. 



The problem of conductivity of the random composite medium and the equivalent problem 
of diffusion in a symmetric (self-adjoint) random environment has been a subject of intensive 
study for the last 25 years. It is virtually impossible to give a full reference list and we just 
mention few papers where the mathematical aspects of the theory were considered for the 
first time: [12], [15], [16], [1]. In the mathematical literature this problem usually is quoted 
as the problem of homogenization for the second order elliptic differential operators with 
random coefficients. Roughly speaking the main result can be formulated in the following 
way: there exists a non-random effective conductivity tensor or effective diffusion matrix 
such that the asymptotic properties of the system are the same as for a homogeneous system 
governed by the effective parameters. The subject is a very active research area till now with 
a vast number of papers publishing every year. However there are very few results related 
to the problem of calculation of effective conductivity and diffusion matrix. In addition to 
the trivial one-dimensional case such results are known only in the self-dual situation in 
dimension two (Kcllcr-Dykhne duality) and in the case of two-component systems where 
the analytic continuation method is used to express the effective conductivity as an analytic 
function of the ratio of the conductivities of two components (see [2], [14], [8], [4], [5]). In this 
paper we discuss a very general rigorous method in the lattice case which was developed in 
[1]. The method is based on a convergent power series expansion for the effective parameters 
and can be applied for arbitrary probability distribution of random conductivities. However, 
the combinatorics of this expansion is rather complicated. That is a reason why it was not 
used for concrete calculations in the past. 

The present paper has two main goals. First of all we demonstrate the constructive 
potential of the method in [1] and give exact formulae for the first 5 orders of the expan- 
sion for the effective conductivity in arbitrary dimension. In the 2D case we also calculate 
the 6th order terms. We then use our exact results to study the quality of the classical 
Bruggeman approximation. We show that in the 2D case the Bruggeman approximation is 
extremely accurate and coincides with the exact answer up to the terms of the 4th order. 
We assume everywhere that the random conductivities (jump rates) are independent iden- 
tically distributed random variables. Although we consider only the case of lattice we 
strongly believe that the method can be generalized for other types of lattices and even for 
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the continuous situation. 

Yakov Sinai was a teacher of one of us and it is our pleasure to dedicate this paper to his 

65th birthday. In fact one of the motivations for this paper was to illuminate the method 
developed together with Yakov Grigorevich and to demonstrate its effective power. 

2 Effective conductivity on Lattice. 
2.1 Exact expansion for effective conductivity. 

We consider effective conductivity for a random wire problem on the cZ-dimensional cubic 
lattice if- ^ d > 2. Throughout the paper we assume that bond conductivities cr are indepen- 
dent identically distributed positive random variables. We are not making any assumptions 
on a probability distribution of a which can be either discrete or continuous. As we have 
mentioned above the calculation of the effective conductivity is equivalent to the calculation 
of the effective diffusion matrix for the continuous time random walk in random environ- 
ment. In this case random conductivities should be understood as jump rates through the 
corresponding bond. We shall use the formula for the effective diffusion matrix which 
was obtained in [1]. This formula is given by a convergent series where the role of small 
parameter is played by a deviation of a random variable a from its average value (a) . Since 
we consider transitions only along the bonds of Z*^ lattice with i.i.d. transition rates a, the 
effective diffusion matrix is a scalar matrix: Me = 2ael, where effective diffusion coefficient 
(or effective conductivity) CTg can be expressed in terms of a convergent power series. We 
first introduce the necessary notations. 

A path 7 = {(-^1, CKi), {z2, 0^2), ■ ■ ■ , {zk, ctfe)} is a finite sequence of pairs {z, a) where 2; is a 
point of lattice Z'^ and a — 1,2, ... ,d corresponds to one of the d possible directions. Notice 
that Zi, Zi+i are not necessarily neighbours on the lattice. The sum of two paths 7 = 71 + 72 is 
simply the ordered union of two sequences where the pairs of the second path follow the pairs 
of the first one. With each pair {z, a) we associate a random variable aa{z) — a{z, z + e^), 
where is a unit vector in the direction a and a{z, z + ea) is the random transition rate 
(conductivity) along the bond {z,z + e^). Denote by Ua{z) = and define for each 
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path 7 = {{zi, CKi), {z2, 0:2), ... , {zk, CKfe)} the moment 

{l)-(tluaM). (1) 

i=l 

A convergent expansion below for the effective conductivity is expressed through the mo- 
ments of a random variable u. We shall also need the following cumulant of a path 7: 

fc m 

^(7) = E(-ir"' E IK^^)' (2) 

m=l 71 H h7m=7j=l 

where summation in (2) is taken over all possible partitions of the path 7 into a sum of paths 
■jj. Finally we define a kernel Tajslz): 

/-^ sin7rAaSin7rA;3COs27r((A,z) - + |A^) -[^ 
ra/3U) = -/ .../ T^- [[dX^, (3) 

where A = (Ai, . . . , A^). Notice that rQ,o,(0) = — | and Tapiz) — r^Q,(— 2;). We can now write 
the following exact formula for ag'- 



<^e = M(l + 5]A«), (4) 

where 



fc=2 



fc-1 



7={(^i-«i)--,fe,«fc)}G£;p^ ^=1 
Here ^j'"'' is the set of all possible paths 7 = {{zi,ai), . . . , {zk.ak)} such that zi = and 
0:1 = 0;^= 1. It has been proven in [1] that the infinite sum in (5) is absolutely convergent. 
That is due to the fact that for the paths 7 which might lead to divergence of A^*^) one has 
£'(7) =0. It was also shown that the expansion in (4) is absolutely convergent and gives an 
exact value of cxe provided \u\ < uq < 1/2. The last condition is technical and probably can 
be improved. In the following proposition we rewrite (4), (5) in a slightly different way. 

Proposition 1 ([1]). 

Assume that there exists a constant mq < | such that \u\ < uo with probability 1. Then for 
any dimension d 



^ 00 [2] 

E <....M')---^^'-)\^ (6) 



fc=2 m=l si,...,sm>2 
siH |-sm=fc 



where the constants aif,,..,sm depend only on dimension d and [•] denotes the integer part. 
Moreover, for any n>l the following estimate holds 

( „. \%\ \ 



I k=2 m=l si,...,sm>2 
\ siH |-sm=fc 

Note that the series in (6) is absolutely convergent. 



< 



(2»o)"+^ 
1 - 2uo 



(7) 



2.2 The 4th order expansion. 

It is easy to see that only those paths for which each pair {z, a) is present at least twice give 
nonzero contribution to (5). This immediately implies that 



d 



Hence the 3-rd order approximation to CTg is given by 



d^ 



(8) 



ar = {a) 1 



(9) 



d d^ 

In the 4-th order the combinatorics is slightly more complicated. Indeed, nonzero contribu- 
tions correspond to the paths 

7(4) = {(0, 1), (0, 1), (0, 1), (0, 1)} , 7t(4) = {(0, 1), {z, 1), {z, 1), (0, 1)}, z ^ 

7a,.(4) = {(0, 1), {z, a), {z, a), (0, 1)}, a^l, 7i%(4) = {(0, 1), {z, 1), (0, 1), {z, 1)}, z ^ 

Another possible type of paths 7i 2(4) = {(0, 1), (0, 1), {z, 1), {z, 1)}, z ^ gives zero contri- 
bution since E{'yf^{4:)) = 0. Easy calculation gives 



+ 



{urrn{o){j2TUz)-rlM 



+ 



a=2 



Notice that 



Hence 



Er?„W=/ •••/ 

^r^^ Jo Jo 



^ '"^ sin^TrAflSin^TrAr 



(E-y=iSin 7rA^)%=i 



Y[dx^ 



E Er?J-) = i- 

/3,a=l z^lA 



(10) 



(12) 
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E E ^U^) (13) 



a=l 



a=l 



Since 



does not depend on (5 we get 



Using (10, 14) we obtain 

^(^) = - '^-^{uY + (urE^Uz) . (15) 

The third term in (15) vanishes in the 2D case. Indeed, ii z — {x,y) we have Ti-i{x,y) — 
^22(y,x). Obviously ru{y,x) + r22{y,x) = if {y,x) ^ (0,0). Hence, for nonzero (x,y) 
we have rii{y,x) = —rii{x,y) which immediately implies Ylz^o^iii^) ~ ^- a result we 
obtain the 4-th order approximation for d = 2: 

= + • (16) 

We next demonstrate that for d > 3 

Ern(-)v^-^ (17) 

which implies 

Ern(^)^0- (18) 
Denote H{d) — —d^Y^zeZ'' l'ii('^)- Using simple Fourier analysis we have 

H(d) = f ... f H{X, /.) JT dX, JT dii, , (19) 

where 

urx .,^ _ sin^(7r(Ai + /xi)) sin^(7rAi) sm^jn^ii) 

dE7=iSin (7r(A^ + /i^)) ^E7=iSin (ttA^) ^ E7=i sm (vr/x^) 

As we have explained above the symmetry in the 2D case gives Ez^ol'ii(^) ~ ^ which is 
equivalent to H{2) — 1. We conjecture that H{d) is a strictly decreasing function of d. The 
conjecture implies that J^z^o^iii'^) ^ ^ d > 3. Although the conjecture above was 

not proven rigorously we have checked it numerically for 3 < d < 5: 

H(3) = 0.923 , i/(4) = 0.874 , H(5) = 0.846 . (21) 



Finally, we get the following 4-th order approximation in an arbitrary dimension: 

r(4) = lrT\(^ - l/,,2\ , _ }_l,A\ _ d+H{d)-2, 



2.3 The 5th order expansion. 

We proceed with the 5-th order calculations. The following paths give nonzero contributions: 

7(5) = {(0, 1), (0, 1), (0, 1), (0, 1), (0, 1)} , 7^^,(5) = {(0, 1), {z, a), {z, a), {z, a), (0, 1)} 
7^,,(5) = {(0, 1), (0, 1), {z, a), {z, a), (0, 1)} , 7^.(5) = {(0, 1), {z, a), (0, 1), {z, «), (0, 1)} 
7a,.(5) = {(0, 1), {z, a), {z, a), (0, 1), (0, 1)} , 7^(5) = {(0, 1), {z, 1), (0, 1), {z, 1), {z, 1)} 
7i%(5) = {(0, 1), (z, 1), {z, 1), (0, 1), {z, 1)} , 7t(5) = {(0, 1), (z, 1), (0, 1), (0, 1), {z, 1)} 
7t(5) = {(0,l),(0,l),(^,l),(0,l),(^,l)}. 

Notice that in the case a — 1 the summation in the paths 7i,2(5), 7f 2.(5), 1 < s < 4 is 
performed over all z ^ 0. Using (5) we get 

A^^) = l-^{u')+K,{d){u'){u'), (23) 

where 

^'»(<') = E E rL(-) + E rtJ^)) - ^ - 5 E ■ (2") 

a=l \ zeZd zeZd / z^O 

This together with (14) gives 

In the 2D case both the first and the last term in (25) vanish and 

^5(2) = Yl rii(^) + E ri2(^) = A + (26) 

where 



Jo Jo 



sin^7r(Ai -/Xi)sinV/iici/ii(i/X2 

/iiiAi,A2j - 7 w Y 

' sin^ 7r(Ai — Hi) + sin^ 7r(A2 — /X2) ) I sin^ tt/^i -|- sin^ 7r//2 ) 
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and 

»i pi 



h ^ f f hl{X,,X2)dX,dX2, 

Jo Jo 

, . . sin7r(Ai - /^i) sin7r(A2 - /X2) sinTT/^i sin7r/i2 (i/ii(iA*2 ^r,o\ 
ri2(.Xi,X2) ^ "7 w \--\^^) 

° ° I sin^ 7r(Ai — /Xi) + sin^ 7r(A2 — A*2) ) I sin^ tt/Xi + sin^ 7r/X2 ) 



The values of /i,/2 were found numerically: Ji = 0.06391, /2 = 0.00439. As a result we get 
in the 2D case the following 5-th order expansion: 

= (^) (1 - li^') + - liu') + ^{u') + I{u'){u')^ , (29) 

where I = + = 0.0683. 

In the general case d > 3 we have 

d d 

E E ria(^) = E rii(^) + E E rL(^) = Hd) + {d- (30) 



where 



h{d) = f ... j\\{X)\{dX,, 
Jo Jo 1 



^ _ /"^ /"^ sin^ 7r(Ai - fii) sin^ tt^i [1^=1 ^3^^ 



^0 



E^=i sin2(7r(A^ - /x^))^ ^^^^ sin2(7r/x^)^ 



and 



1 „1 d 



I,(d) = f ... f hl{X)fldX^, 
Jo Jo 



Jo 



/ /-i sin7r(Ai -/ii)sin7r(A2 -/i2)sin7r/iisin7r;U2 
/12(A) = / . . . / 7 w 4 ■ (32) 



E^=i sin2(7r(A^ - /x^))^ ^^^^ sin2(7r//^)^ 



Collecting all the terms we get 



where I{d) = /i(ci) + (d - l)l2{d) and //(d) is given by (19), (20). 



2.4 Keller-Dykhne duality and the 6th order expansion in the 2D 
case. 

Although it is possible in principle to calculate an expansion of an arbitrary order the problem 
becomes more and more cumbersome for higher order terms. However in the 2D case one 
can significantly simplify calculations using the duality symmetry which was discovered by 
Keller ([9]) and Dykhne ([6]). Consider duality transformation 

a^\. (34) 

Denote by {a}, {cr^^} the probability distributions for positive random variables u and 
respectively. Then duality symmetry which holds only in the 2D case implies that 

a.iW}) = a:\{a}) . (35) 

Although both Keller and Dykhne considered only the continuous systems the symmetry 
(35) can be extended to the case of discrete lattice systems which we study in this paper 
(see [11]). The duality symmetry immediately implies that in the self-dual case. i.e. when 
the probability distributions {cr} and coincide, the effective conductivity — 1. 

It also gives an exact answer in the case which we call almost self-dual. We say that the 
probability distribution for a random variable a is almost self-dual with respect to the duality 
transformation (34) if there exists a positive constant (Tq such that the probability distribution 
for (Jocr is exactly self-dual, i.e. 

{aoa ] = {{aoa)-^]. (36) 

Since Ue is a homogeneous function of the first order and ae({o"o(j}) = 1, it follows that 
in the almost self-dual situation (Te({cr}) = a^^ . Notice that in the two-component case 
with equipartition, i.e. when a takes values Ui and (72 with probabilities | the probability 
distribution for a is almost self-dual with ctq = {y/o\a2)~^. Hence, 

(7e = ctq"^ = ^/cr . (37) 

This well-known result by Keller and Dykhne provides one of the very few exact solutions 
for the effective conductivity. 

We next show that the duahty symmetry alone gives a lot of relations on the coefficients 
of the expansion (6). In fact we shall be able to recover the 6th order expansion using only 
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the 5th order and the symmetry. Consider the case when a takes three values: 1 — e with 
probabihty p, \ — ae with probabihty p and 1 with probabihty 1 — 2p. Correspondingly a 
random variable g^^ takes values and -^^^ with probabilities p and 1 with probability 
1 — 2p. We shall use the formula (6) in order to calculate cre({cr})(Je({cr~^}) and check the 
duality identity (35) subsequently in the 2nd, 4th, 6th and 8th orders of the power series 
expansion in e. This inductive procedure allows to find all the relations on the coefficients 
a>su...,s^- We performed calculations using the Maple symbolic package. In the 2nd order one 
immediately gets a\ = The 4th order calculations give two relations: 



2 

3 (2) ^ 1 
2'-' S'""^ ~4 2 



(2) -5 (2) <3 (2) -L O (2) , , 

«2,2 = o«3 - o ' «4 = 7 - o«3 ■ (38) 



The 6th order expansion provides four more relations: 

7 . 3 (o\ 15 C2\ 1.1 



(2) - LJ^) ^ ^J'^) _ „(2) _ ^^^(2)^2 _o^(2) _^(2) 

3 "2,3 > 



"2,2,2 — r)"3 ^ o 2,3 -.^ 1 "3,3 — r, ^ O ^ 3 J ^" 



20 2 ^'•^ 16 ' 2 2 

(2) 11 ^ (2) 3 (2) 5 (2) (2) 5 (2) 5 (2) 1 , . 

«2,4 = y - 64 ' - 2^2,3 + 2^5 ' "6 = 2'^3 ~ 2^^^ ~ 2 ■ ^ ' 



Using (29) we have 

af = \ , = ^,a?], = I^ 0.0683 (40) 
which immediately gives a'^^ — , af^ — and 

J^) fl(2)_J__. „(2)_J__3. ^(2)__J_ 

«2,2,2 - 2^ 16 ' «3,3 - 32 ^ ' «2,4 - 32 3 ' ^ ~ 32 ' ^ ^ 

As a result we obtain the 6th order expansion in the 2D case: 




3 The Bruggeman Approximation. 
3.1 Bruggeman's equation. 

The Effective Medium Approximation (EMA) was invented by Bruggeman [3], and has re- 
mained one of the most popular approximations used for calculations of the linear bulk 
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effective electrical conductivity (7e of a many-component composite medium. This is mainly 
due to the simphcity of EMA and to the fact that it gives accurate results for a wide range of 
parameters. It also has a non-trivial percolation threshold which most other simple approx- 
imations do not possess. Another advantage of Bruggeman's approximation is connected 
with the fact that none of the complicated details of the microstructure are used in its con- 
struction. EMA is only based on the assumptions that the composite is macroscopically 
homogeneous and isotropic and that individual grains are spherical. It is also important 
to mention that EMA applies without any changes to the calculation of dielectric suscepti- 
bility, magnetic permeability, thermal conductivity and chemical diffusion coefficients, since 
in all those cases the mathematical structure of the equations is the same as for electrical 
conduction. 

Suppose that the values of the component conductivities ctj and the component volume 
fractions pi are given. Then Bruggeman's equation in the ci— dimensional case has the fol- 
lowing form: 

n 

Yp, ^0. (43) 

1=1 ^ ' 

^From the mathematical standpoint it has many beautiful properties which are of high 
importance for the theory of random composites. Equation (43) has a unique positive root 
CB((7i) which is homogeneous of the 1-st order, monotone and reducible with respect to the 
equating of some constituents. It is also 5'^-permutation invariant in the case when all are 
equal and compatible with a trivial solution gb = o when all Oi = a . Finally, in the case 
d — 2 the Bruggeman's solution is self-dual with respect to the duality transformation (34). 
Namely, if (7j a^^ and Pi are unchanged then 

(7s((jf \ a^^, a-^) = cT^^((7i, 0-2, (7„) . (44) 

It follows that aB coincides with Keller-Dykhne solutions in the self-dual and almost self-dual 
situations. In particular, as — \^o\a2 for the two-component system with equipartition and 
conductivities taken values ai,a2- Notice that the Bruggeman approximation is also exact 
in the ID case. 
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3.2 Solution of Bruggeman's Equation. 

Let cr be a random variable corresponding to random conductivity. Then Bruggeman's 
equation (43) can be written in terms of averages in the following form 

^7^^, \=0. (45) 

Notice that (45) is the most general form of Bruggeman's equation. We first show that 
Bruggeman's equation (45) has a unique positive solution ub- Indeed, function 

is obviously decreasing. Also F(0) = 1 and F(x) — > — ^ as x — > oo which implies the 

existence and the uniqueness of the solution. We next find the expansion of CTb in terms 



dimensionless variables 



of the moments of the disorder parameter u = ^-j^- It is convenient to introduce new 



Obviously u = r] — 1. In the new variables Bruggeman's equation (45) takes the form 

~ ^ \ - , (4J 



where 6 — d — 1. Notice that 



After the averaging of the both sides in (49) we get 



which together with (m) = immediately implies 

If the random variable u is small enough the solution of equation (51) can be written as a 
convergent expansion in terms of the moments of u: 

oo [i] 

^ = 1 + EE E bf)_sM^)...{u^-). (52) 

k=2 m=l si>--,sm>2 
siH |-sm=fc 
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Notice that this expansion has similar structure to the expansion (6). Easy calculation leads 
to the following expansion up to the terms of 6th order: 

+ ^-^{u'){u') - - ^(«')^') (53) 

2c?- 3, 3 2 2d'^-8d+7 



which gives the 6th order approximation for the Bruggeman approximation 

= (54) 

and its 2D version 

of = (o) - \(^') + j<»'> - \(-^') + ]^<»=> + ^(»^)<»=) 

-^(^V^{«^>(»V^(»V + ^("V). (55) 

3.3 Effective conductivity and the Bruggeman approximation. 

It follows from (33), (53), (54) that the Bruggeman approximation ub coincides with the 
effective conductivity Ue up to the terms of 3rd order. However if d > 3 the 4th order terms 
are different. Let us assume that 

|«| < e , {u^) > ce^ . (56) 



Then, 



'ye- OB- (a) C-^{uY + 0{e')) > {a) ('-^ch^ + 0{e>) ] . (57) 



This implies that for e small enough (Te> (Tb- In the 2D case the Bruggeman approximation 
is even more accurate. It coincides with a^. up to the 4th order terms. Nevertheless, if (cr^) 
does not vanish then Ue differs from (Tb in the 5th order. Assume that (56) holds and in 
addition \{u^)\ > ce^ ■ Then, 

a,-aB = (a) Ul - ^) {u'){u') + 0{e')) . (58) 
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Since / > ^ we have Ue 7^ (Tb for e small enough. Notice that ub is bigger than if {v?) is 
negative. Finally we consider the symmetric 2D case. We shall assume that u satisfies (56) 
and {u^) = 0. Then, 



It follows from (59) that (7e < (Jb if {{u^ — {u^))"^) is of the order of and e is small enough. 
We summarise all three cases in the following simple proposition. 

Proposition 2 

1. Consider the case d > 3. If u satisfies (56) then there exists e{d,c) > such that 
(Je > (Tb for all e < e{d, c). 

2. Let d = 2, u satisfies (56) and \{u^)\ > ce^ . Then there exists e(c) > such that 
7^ CTfi for all e < e{d, c) and sgn (a"e — (Jb) = sgn {{u^))- 

3. Let d — 2 and (u^) —O.Ifu satisfies (56) and {{u'^ — (u^))^) > ce^ then there exists 
e(c) > such that < gb for all e < e(c). 

The following corollary follows easily from Proposition 2. Consider the n— component system 
where a takes the values ui, (J2, . . . , cxn with probabilities pi,P2, • • • ,Pn, Pi > 0, Pi + P2 + 

■ ■ ■ + Pn = 1- Wc shall also assume that the system is irreducible, i.e. cXj 7^ aj, 1 < i, j < n. 
Denote Pr^m = min(pi,p2, • • • ,Pn)- 

Corollary 1 

1. Let d> 3. If \u\ < e((i,pmin) then > ub- 

2. Let d — 2. Assume that n — 2 and pi > p2. Then 0^. ^ ub provided \u\ < e{c), where 
c — P2 (^^ — )• Moreover, if ^2 > (Ti then > (Tb- In the opposite case, i.e. when 
(7i > (72 one has ctb > Ce- 

3. Let d = 2. Assume that n = 3 and (u^) = 0. Then there exists Ci{pi,p2,P3) > 
such that (Te < (Tb if \u\ < e(ci). In particular, if ai = 1 + e, (T2 = IjC^s = 1 — e and 
Pi — Ps — P, P2 — ^ — "^P, < P < \ then Ci(pi,p2,P3) = 2p(l — 2p) and < (Jb under 
condition \u\ < e(2p(l — 2p)). 

Finally, we conjecture that for n— component systems the efi^ective conductivity coincides 
with the Bruggeman approximation only if the probability distribution {cr} is almost self- 
dual, see (36). 




(59) 
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4 Concluding remarks. 



1. We have derived the exact formulae for the first 5 orders of the expansion of the effective 
conductivity in terms of the moments of the disorder parameter u in arbitrary dimension. 
In the 2D case we have also found the 6th order terms. It is quite interesting to extend these 
results to other types of 2D lattices and to the continuous plaquctcs systems. Notice that 
our duality analysis holds for the general 2D case. Hence, if the expansion (6) is valid, it is 

(2) (2) (2) 

enough to find 03 , ag , 03 3 in order to determine all other terms up to the 6th order. 

2. We have shown that Bruggeman's solution (55) gives a remarkably accurate approximation 
for the effective conductivity of the 2D random many-component lattice wire system. It turns 
out that in the case of square lattice the first four orders of the expansion of Bruggeman's 
solution in terms of the moments of the disorder parameter coincide with the corresponding 
expansion of the exact solution. However, in the 5th order the Bruggeman approximation 
deviates from the exact one. An interesting and natural problem is to verify whether such 
behaviour is characteristic for the square lattice or it also holds for other 2D lattices. It is also 
interesting to analyse the relation between Bruggeman's solution and effective conductivity 
for the continuous 2D random composites. Recently four isotropic three-component S^- 
permutation invariant regular structures with three-fold rotation lattice symmetries in the 
2D case were treated numerically [7] . A simple cubic equation with one free parameter A> 

3 

(7^ + A JiCTg - A J2(7e - J3 = , = ^ (Jj , J2 = ^ (7j(7j- , J3 = (7i (72(73 

i=l iytj 

was proposed as an algebraic equation of minimal order. Its solution share many properties 
with (Te and corresponds to Bruggeman's solution when A = ^. The numerically estimated 
values of A corresponding to different cases were calculated with a very high precision. 
It appears that they are distinct and lie rather far from | for some of the structures. This 
indicates a strong dependence of (jg on plane symmetries in contrast with the two-component 
case. 

3. Recently in the paper by A. Kamenshchik and I. Khalatnikov ([10]) the perturbation the- 
ory was developed for the periodic three-component plaquetes lattice systems with two-fold 
rotation lattice symmetry. We hope that their technique combined with our approach will 
lead to the exact expansion for the effective conductivity in the random plaquetes situation. 

4. After the paper was submitted we were informed about the paper by Jean-Marc Luck 

15 



([13]) where very similar results were obtained using different method for calculating an 
expansion for the effective conductivity. In our opinion the approach we use has certain 

advantages. First of all, it is rigorous and, hence, more suitable for mathematical audience. 
Secondly, it gives arbitrary good rigorous bounds for the effective conductivity (see (7)). 

5 Acknowledgment. 

The authors would like to thank A. Kamenshchik and I. Khalatnikov for useful discussions. 
We are also grateful to Jean-Marc Luck for bringing his paper ([13]) to our attention and to 
D. Khmelev for his help with numerical calculations. 

The main part of the paper has been written during the stay of one of the authors (LGF) 
at the Isaac Newton Institute for Mathematical Sciences and its hospitality is highly appre- 
ciated. 

The work was supported in part by grants from the Tel Aviv University Research Authority 
and the Gileadi Fellowship program. 

References 

[1] V. V. Anshelevich, K. M. Khanin and Ya. G. Sinai, Symmetric random walks in random 
environments. Comm. Math. Phys. 85: 449 (1982). 

[2] D. J. Bergman, The dielectric constant of a composite material - a problem in classical 
physics, Phys. Rep. C 43: 377 (1978). 

[3] D. A. G. Bruggeman, Berechnung verschidener physikalischer konstanten von heteroge- 
nen substanzen, Ann. Physik (Leipzig) 24:636 (1935). 

[4] O. Bruno and K. Golden, Interchangeability and bounds for the effective conductivity 
of the square lattice, J. Stat. Phys. 61: 361 (1990). 

[5] A. R. Day and M. F. Thorpe, The spectral function of random resistor networks, J. 
Phys. 8: 4389 (1996). 

[6] A. M. Dykhne, Conductivity of a two-dimensional two-phase system, Sov. Phys. JETP 
32: 63 (1970). 

16 



[7] L. G. Fel, V. Sh. Machavariani and D. J. Bergman, Isotropic conductivity of two- 
dimensional three-component symmetric composites, J. Phys. A 33:6669 (2000). 

[8] K. Golden and G. Papanicolaou, Bounds for effective parameters of heterogeneous media 
by analytic continuation. Comm. Math. Phys. 90: 473 (1983). 

[9] J. B. Keller, A theorem of the conductivity of a composite medium, J. Math. Phys. 5: 
548 (1964). 

[10] I. M. Khalatnikov and A. Kamenshchik, A diagram technique for perturbation theory 
calculations of the effective conductivity of two-dimensional systems, JETP 91: 1261 
(2000). 

[11] P. Kogut and J. Straley, The bicritical macroscopic conductivity exponent in dimen- 
sionalities two, three and infinity, J. Phys.C 12:1 (1979). 

[12] S. M. Kozlov, The averaging of random operators, Mat.Sb. 109(151): 188 (1979). 

[13] J. M. Luck, Conductivity of random resistor networks: an investigation of the accuracy 
of the effective-medium approximation, Phys. Rev. B 43:3933 (1991). 

[14] G. W. Milton, Bounds on the transport and optical properties of a two-component 
composite material, J. Appl. Phys. 52: 5294 (1981). 

[15] G. Papanicolaou and S. R. S. Varadhan, Boundary value problems with rapidly oscillat- 
ing coefficients. In: CoUoquia mathematica societatis Janos Bolyai 27, Random fields, 
V.2, North-Holland Publ. (1981). 

[16] V. V. Zhikov, S. M. Kozlov and O. A. Oleinik, G-convergence of parabohc operators. 
Russian Math. Surveys 36:11, (1981). 



17 



